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Can X-rays provide a solution to the abundance discrepancy 
problem in photoionised nebulae? 
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ABSTRACT 

We re-examine the well-known discrepancy between ionic abundances determined via the 
analysis of recombination lines (RLs) and coUisionally excited lines (CELs). We show that 
abundance variations can be mimicked in a chemically homogeneous medium by the presence 
of dense (hh ^ 10"* cm"^) X-ray irradiated regions which present different ionisation and 
temperature structures from those of the more diffuse medium they are embedded in, which is 
predominantly ionised by extreme-ultraviolet radiation. The presence of X-ray ionised dense 
clumps or filaments also naturally explains the lower temperatures often measured from O 11 
recombination lines and from the Balmer jump when compared to temperatures determined 
by CELs. We discuss the implications for abundances determined via the analysis of CELs 
and RLs and provide a simple analytical procedure to obtain upwards corrections for CEL- 
determined abundance. While we show that the abundance discrepancy factor (ADF) and the 
Balmer Jump temperature determined from observations of the Orion Nebula can simulta- 
neously be reproduced by this model (implying upward corrections for CELs by a factor of 
1.15), we find that the required X-ray fluxes exceed the known Orion's stellar and diffuse 
X-ray budget, if we assume that the clumps are located at the edge of the blister We propose, 
however, that spatially resolved observations may be used to empirically test the model, and 
we outUne how the framework developed in this letter may be applied in the future to objects 
with better constrained geometries (e.g. planetary nebulae). 

Key words: galaxies: abundances, (ISM:) H ii regions, (ISM:) planetary nebulae: general, 
ISM: abundances. X-rays: general 
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1 INTRODUCTION 

The analysis of emission line spectra from photoionised gas (e.g. 
in H II region and planetary nebulae, PNe) often provides one of 
the best method to determine elemental abundance in our own and 
other galaxies. Accurate metallicities are fundamental to the char- 
acterisation of many astrophysical system, and a number of em- 
pirical techniques have been developed in order to extract this in- 
formation from the observation of nebular emission lines. Histori- 
cally, coUisionally excited lines (CELs) provided the best options, 
thanks to their brightness and their temperature and density diag- 
nostic power. One of the major drawbacks of CELs is their expo- 
nential dependence on electron temperatures, which introduces un- 
certainties in the derived abundances, particularly for high metal- 
licity objects, where large temperature gradients may occur (e.g. 
Stasinska et al. 2005). 

Recombination lines (RLs), produced by electron cascades 
following recombination of metal ions, are potentially more reli- 
able abundance estimators, since they have a much weaker (inverse) 
dependence on temperature and are also only weakly dependent 
on density variations in the nebular regime. However, due to their 
faintness (roughly 10"^^ relative to Hfi), they are only detectable in 



nearby objects and through very deep observations. Nevertheless 
RLs from e.g. O II and C II are now routinely detected in many 
bright H II regions in the Galaxy and in the Magellanic clouds 
(e.g. Tsamis et al. 2003; Garcia-Rojas & Esteban, 2007, and refer- 
ences therein) , as well as in many PNe (Tsamis et al. 2004, 2008; 
Garcia-Rojas et al. 2009; Wesson et al. 2005, 2008; Wang & Liu, 
2007; Krabbe & Copetti, 2006; Ercolano et al. 2003b, 2004). How- 
ever, already over sixty years ago, comparison of ionic abundances 
obtained by RLs and CELs had revealed a systematic discrepancy 
such that RL-determined abundances are always found to be higher 
than those determined from CELs (Wyse 1942; Aller & Menzel, 
1945), these findings have been confirmed for all observed nebu- 
lae. The problem affects both H II region abundances, with abun- 
dance discrepancy factors (ADF) of order two (e.g. Garcia-Rojas et 
al. 2007, and references therein), and PNe, where very large ADFs 
have been reported for some extreme objects (e.g ADF -80 for 
Hf2-2, Liu et al. 2006). This is known as the abundance discrep- 
ancy problem and is probably the major outstanding problem of 
nebular astrophysics (for further discussion see Sections 5.11 and 
9.5 of Osterbrock & Ferland, 2006). 

The most important question is which set of abundances is 
more reliable, a question that can only be answered when the cause 
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of the discrepancy is known. A number of solutions have been pro- 
posed, which include (i) temperature fluctuations, also known as 
the paradigm (e.g Peimbert 1967), and (ii) hydrogen deficient in- 
clusions (e.g. Liu et al. 2000; Stasihska et al. 2007). It is beyond 
the scope of this letter to review the extensive research carried out 
by the supporters of (i) or (ii), however it should be noted here 
that temperature fluctuations in chemically homogeneous objects 
imply that CELs underestimate abundances and that RLs provide 
the more accurate estimate. The presence of metal-rich (hydrogen 
deficient) inclusions, implies instead that RLs are predominantly 
emitted from these inclusions and therefore the abundance derived 
from RLs is strongly affected by this gas-phase and it is not repre- 
sentative of the bulk of the nebula. 

Both theories have limitations, and a large body of literature 
exist in support and against one or the other. Problems with the t' 
paradigm include the fact that it cannot be used to explain ADFs 
larger that a few, furthermore photoionisation models have so far 
failed to reproduce the temperature fluctuations required to explain 
even low values ADFs (e.g. Ercolano, Bastian & Stasihska, 2007). 
A problem with the metal-rich inclusions is in proving their exis- 
tence, a number of origins have been proposed for H II regions (e.g. 
Stasihska et al. 2007) and PNe (e.g. Ercolano et al. 2004), but no 
conclusive evidence has been brought forward. 

In this letter we investigate whether the solution to the abun- 
dance discrepancy problem may lie in the existence of X-ray irradi- 
ated high density (> lO'^cm"^) quasi-neutral clumps and filaments 
embedded in the more tenuous EUV-ionised nebular gas. We show 
how this chemically homogeneous two-phase medium may natu- 
rally account for typical H II region ADFs, and discuss existing ob- 
servational evidence in support of our theory. We conclude however 
that the known X-ray (diffuse and point) sources in H II regions, fall 
short of the requirements of our model. 

In Section 2 we derive a relation of the ADFs to the rela- 
tive electron temperature and ionisation structure of the two phases 
and their relative ionised masses. In Section 3 we discuss the pho- 
toionisation process in a quasi-neutral gas irradiated by X-rays and 
present results of photoionisation modelling of plane parallel slabs 
using typical conditions, as inferred from available observations. In 
Section 4 we further discuss our results and their implications for 
the determination of gas abundances from nebular emission lines in 
the Orion Nebula. Our conclusions are summarised in Section 5. 
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where the E and X superscripts denote physical quantities for the 
EUV and X-ray ionised regions (regions E and X), respectively, 
and we have defined the ionised hydrogen emission measure as 



EM = x(//+) ■ iiH -n.-V 



(2) 



We note that, since we are considering a chemically homogeneous 
medium, A/H is the same in both regions. 

We define the constant e which represents the ratio of the 
ionised hydrogen emission measures in the two regions: 

= E ■ EM^ (3) 

RL emissivities carry an approximately inverse dependence 
on electron temperature and are also weak functions of electron 
density, therefore we define i and t' such that. 
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Substituting Equations [3] to |4] into Equation [T] and rearranging we 
obtain 
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where we have defined 7 = ^ - . 

is what an observer would use to determine ( 4 ■4§t^ I from 
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H 'mW) ) ' where the CEL analysis is based on the flux ratio 

of CEL lines (e.g. [O///]/l5007, 4959 A) to Hji at a given T^. As 
shown in the next two sections, temperatures in a quasi-neutral X- 
ray heated region (X region) are lower than those in the ionised 
EUV heated region (E region) and therefore we can assume that the 
CEL flux is dominated by region E (due to the exponential depen- 
dence on electron temperature of the CEL emissivities). However 
Hp is likely to be produced in both regions, in proportions dictated 
by the parameters i and e, defined in Equation[3] Namely, 



2 ABUNDANCE DISCREPANCY FACTOR IN A 
TWO-PHASE MEDIUM 

We consider the case of a two phase medium comprising high 
density (>10'' cm"'') clumps ionised by X-ray radiation embed- 
ded in a lower density component (e.g. <1000 cm"''), predom- 
inantly ionised by the extreme-ultraviolet (EUV) radiation field 
from young OB stars, clusters or complexes. In the next section 
we will briefly review photoionisation processes in quasi-neutral 
X-ray irradiated gas, but for the sake of the current discussion 
it is sufficient to define the ion fraction of ion A' of species A 
as x(A'), such that the total number of A' ions in a region with 
neutral hydrogen number density and volume V is given by 
N(A') = x(A') ■ {AIH) ■ tiH ■ V, where A/H is the total abundance 
of element A by number relative to hydrogen. We further define 
juL as the temperature and density dependent emissivity of a given 
emission line produced by recombination of A' + e" — » A'"' + hv. 

The observed RL/Hp flux ratio (%rl) emerging from such a 
region is given by 



I(Hp)"" = iJlp ■ EM") + (ill, ■ EM") 
= f^^.EM"-{^). 

It follows that 
A x{A') y*' _ A x(A')^ e 

H x(H+) j„, ~ HIqFy 



(6) 



(7) 



and 



(. \obs I . \( 



RL 



obs 
CEL 



where, for conciseness, we have defined the parameter 
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The upward correction for CEL-determined A' abundances, 
;f(A'), can be obtained from equation [T] 
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3 lONISATION STRUCTURE OF DENSE GAS 
IRRADIATED BY A TENUOUS X-RAY FIELD 

Multiwavelength studies of star-forming regions like Orion and 
30 Dor have revealed a highly fragmented structure where warm 
molecular gas and ionised gas are intermixed (Poglitsch et al. 
1995). The high density clumps and filaments are opaque to EUV 
photons, but they are maintained at a low level of ionisation by X- 
rays, which can penetrate much larger neutral hydrogen columns 
(e.g. lO^^cm"- at ~1 keV). The ionisation structure of these re- 
gions is very different from that of a typical H II gas. X-rays can 
photoionise the inner shells of atoms and ions, from which Auger 
multi-electron ejection may follow. This process couples non adja- 
cent stages of ionisation. Furthermore high energy radiation may 
also eject suprathermal electrons which then produce secondary 
ionisations (e.g. Xu & McCray, 1991). 

We have used the photoionisation code MOCASSIN (Ercolano 
et al. 2003, 2005, 2008) to calculate the ionisation and temperature 
structure of dense slabs irradiated by a thermal X-ray spectrum. 
Table 1 lists the model parameters for the illustrative cases pre- 
sented here. The slabs are irradiated using a synthetic spectrum ob- 
tained for an X-ray temperature of log(rx) = 6.8 K using the proce- 
dure detailed in Ercolano et al. (2008). We use "solar" abundances 
from Grevesse & Sauval (1998) with the exception of C,N and O 
abundances which are taken from Allende-Prieto et al. (2002), Hol- 
weger (2001) and Allende-Prieto, Lambert & Asplund (2001), re- 
spectively. The density of clumps and filaments in real nebulae is 
likely to vary, however their density must be large enough to be- 
come opaque to EUV radiation. Poglitsch et al. (1995) suggest den- 
sities as high as 10* cm"'', in agreement with observations of molec- 
ular gas in ionised regions by other authors who inferred densities 
rangmg from 10*'cm-^ to 10* cm . Therefore here we experiment 
with slabs of densities in this range, adjusting the length of the slabs 
such that they are truncated at a column of lO-'cm"^. 

For the range of densities and X-ray fluxes covered by our 
models, we find that very low electron temperatures can be pro- 
duced in the slabs while a low level of ionisation is still maintained. 
Table 1 lists the values of < T,{H*) >, < x{H*) > and < x{0-*) >, 
defined as 



< TAA') >-- 



< x(A') >= 
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The X-ray irradiated slabs show very different characteristics from 
those of EUV irradiated H II gas, which has generally much higher 
electron temperatures (roughly 7000 to 12000 K) and < x(0^+) > 
to < x{H*) > ratios always less than unity (Stasihska, 1982). 

However before proceeding to our discussion of the impli- 
cation of our models for the ADF in Orion, it is worth men- 
tioning a few caveats, (i) The results presented here are for low- 
ionisation/neutral gas and do not include a chemical network for 
reactions typical of X-ray dominated regions (XDR; e.g. see Mei- 
jerink & Spaans, 2005). As a consequence heating and cooling 
channels relevant to these regions have also been ignored. However 
our models are restricted to low column densities (lO^'cm^^) where 
molecular hydrogen is unlikely to be present, (ii) Due to the large 
neutral hydrogen and helium abundances, charge exchange plays 



an important role in the ionisation balance of heavy elements. The 
accuracy of the ionisation structure predicted by our models finally 
depends on the uncertainties in the charge exchange rates adopted 
(Ferland et al. 1997). (iii) The dense clumps are also bathed in a far- 
ultraviolet radiation (EUV) field. While this will not affect the ioni- 
sation structure of the gas, it may affect its temperature, principally 
via photoelectric emission from the surface of grains, (iv) Finally, 
grain depletion may also be more important in dense clumps and 
in the photon dominated region, modifying ADEs for carbon and 
oxygen. In these cases there may be variations within an H II re- 
gion depending on how much of each phase the line of sight passes 
through. The blister geometry of Orion may be a good example. 
These processes are not included in the current models. 

While, in view of these limitations, we caution that the tem- 
peratures reported here are not to be considered accurate, we stress 
however that the 'real' temperatures will still be factors of several 
below the EUV-heated gas temperatures, and therefore the general 
conclusions of our work remain unaffected. Nonetheless this un- 
derlines the need for further modeling which takes into account of 
XDR physics and EUV irradiation. 



4 APPLICATION TO THE ORION NEBULA 

Abundance discrepancy factors (ADEs) have been measured in a 
number of H II regions and cluster around the value of two (e.g. 
Garcia-Rojas & Esteban, 2007 and references therein ). Due to 
its proximity, Orion is one of the best studied regions and has a 
reported mean ADF value of 1.4 (e.g. Garcia-Rojas et al. 2007). 
Mesa-Delgado et al. (2008) presented deep long slit spectra at sev- 
eral positions in the Orion Nebula. They found that ADF{0^*) was 
roughly constant along the slit positions, except at the location of 
Herbig-Haro objects where the largest increases were recorded. In 
the framework of the clumpy X-ray + diffuse EUV (X-E) model, 
the enhanced ADF{0^*) at these locations is due to X-rays from 
the Herbig-Haro shock front impinging on cold neutral material 
(Pravdo et al. 2001). Also in agreement with this model, enhanced 
values of ADF{0^*) were also reported by these authors in the re- 
gion closer to 0' Ori C (which emits ~60% of the entire X-ray flux 
in the Trapezium). 

The physical conditions that can create a mean ADF value of 
1.4 in Orion are now explored in the context of the X-E model. 



We use ( = 



=< r^(0^+) > and IV = 



- /- 



and assume values for the EUV ionised region as determined by 
observations (e.g. Garcia-Rojas et al. 2007) and detailed photoion- 



isation modeling (Baldwin et al. 1991) which give 



-v(0"+)^ 



= 0.57 



and rf = 8340 ± 130 K. The emissivities of the On VI multiplet 
are calculated using the data of Storey (1994) and Smits (1991). We 
use the above to determine t and W from our models (see Table 1) 
and using these in Equation [S] we find the values of e, the ratio of 
ionised hydrogen emission measures in regions E and X, needed to 
produce an ADF of 1.4. From the values of e and l thus calculated, 
X{0^*), the upward correction for the CEL-determined O^* abun- 
dance, can also be easily determined (Equation 11) and is listed in 
Table 1. The O^^ abundance in the E-region, in this case, can be 
obtained from the GEL value after an upward correction by factors 
of 1.12 to 1.27. 

The values of e shown in Table 1, are encouraging as they 
show that the ionised emission measure required for the X region is 
always significantly smaller than that in the E region. An observa- 
tional constraint on e can be teased out from a comparison of tem- 
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Model 


log(Fx) 


log(«H) 






< .v(02+) > 




e'' 


XiCELY 


T(BJ)'' 




[phot/s/cm^^] 


[cm"^] 


[IO'K] 


[10-3] 


[10-3] 








[lO^K] 


F6N4 


6 


4 


1.2 


0.71 


0.90 


6.91 


38.4 


1.18 


5.1 


F6N5 




5 


1.3 


0.22 


0.22 


6.27 


27.5 


1.23 


4.8 


F6N6 


6 


6 


1.2 


0.05 


0.06 


7.24 


26.5 


1.27 


4.2 


F7N4 


7 


4 


1.6 


2.5 


4.1 


5.27 


37.9 


1.14 


6.1 


F7N5 


7 


5 


1.5 


0.83 


0.97 


5.62 


28.8 


1.19 


5.4 


F7N6 


7 


6 


1.4 


0.23 


0.26 


6.02 


29.8 


1.20 


5.1 


F8N4 


8 


4 


2.3 


9.2 


18.2 


3.66 


31.8 


1.12 


6.9 


F8N5 


8 


5 


2.0 


3.3 


4.6 


4.21 


25.7 


1.16 


6.2 


F8N6 


8 


6 


1.8 


0.86 


1.2 


4.68 


28.6 


1.16 


6.0 


F9N4 


9 


4 


4.5 


36 


57 


1.87 


13.0 


1.14 


7.7 


F9N5 


9 


5 


3.2 


12 


20 


2.63 


19.2 


1.14 


7.2 


F9N6 


9 


6 


2.8 


4.1 


6.2 


3.01 


20.0 


1.15 


6.9 



Table 1. Model parameters and results for X-ray irradiated quasi-neutral slabs. " Estimated for Orion assuming '^^u+^e = 0.57 (Balwin et al. 2001), =< 

TeiH*) > and Tf = 8430 K (Garcia-Rojas et al. 2007). The emissivities of the Oil VI multiplet are from Storey (1994) and Smits (1991). *Infen-ed from 
Equation [8] for the Orion ADF=:1.4 (Garcia-Rojas et al. 2007). ^Inferred from Equation 11. ''inferred from Equation 1 131 the observed Tj(BJ) in Orion is 
7900±600 K. 



perature measurements from the Balmer Jump (Tg(BJ)) and from 
coUisionally excited lines (Te(CEL)). Balmer Jump (BJ) measure- 
ments will contain contributions from both the X and the E regions, 
while CELs are expected to be dominated by the E region. Starting 
from this assumption one can derive the following relation: 

r,(Sy)= — ^ -T.iCEL) (13) 

where we have used the fact that the BJ to Hll line ratio is propor- 
tional to temperature to the power of -2.3 (e.g. Liu et al. 2001). 

From the equation above and the values of e and i inferred 
from our models and Orion-type E-region conditions and ADF val- 
ues, we have calculated the predicted Te(BJ) (last column of Ta- 
ble 1), which we can compare with the observed value of 7900±600 
K (Garcia-Rojas et al. 2007). The comparison shows that models 
with an X-ray flux of lO'phot/sec provide the best match, which 
implies an upward correction to the CEL abundance of 1.14. 

4.1 Origins of tlie X-ray flux 

H II regions present a number of sources of X-ray radiation. In the 
Orion Nebula we can simplistically list interacting stellar winds as 
a source of the diffuse component (e.g. Guedel et al. 2008) and 
young stars emitting X-rays as point sources (e.g. Feigelson et al. 
2005). The observed diffuse X-ray source in the Orion is quite 
small (Lx = 5.5 x lO^'erg/s; Guedel et al. 2008) and cannot pro- 
duce fluxes in the range discussed in Section 3. However the stel- 
lar wind itself carries a much larger mechanical luminosity. The 
wind from 6'Ori C alone has L^ech = TxlO^^erg/s. If all of its en- 
ergy were transformed into 300eV photons upon impact with the 
clumps, this would correspond to a photon flux of 1.3x10' phot/s 
for clumps located at the wall of the blister (roughly 0.1 pc dis- 
tance from 6)'Ori C, Baldwin et al. 1991; Schiffer & Mathis, 1974), 
matching the requirements of the X-E model. Considerations of 
radiative timescales, however, make the above rather unlikely. In- 
deed the low density stellar wind will form a bow shock around the 
clump, but, due to the extremely low number density of the wind 
the radiative timescales will be very long and only a very small 
fraction of the impacting mechanical flux will be radiated in the 
X-ray at the shock-front. If the clumps are embedded in the H II 



region gas the wind would have to first interact with the ionised 
material, further complicating matters. 

The ONC stars radiate approximately 5.7xl033erg/sec in the 
X-ray (Feigelson et al. 2005), which corresponds to a photon lumi- 
nosity of roughly 3.5x1 0'^^phot/s (assuming 1 keV photons given 
that the stellar spectra are much harder than shock-produced emis- 
sion). A photon flux of ~10'phot/s, which gives the best fit to 
Orion's ADF and BJ temperature (see Table 1), would require the 
clumps to be placed at a distance of ~1.7xl0'^cm from the source, 
which is about 15 times smaller than the estimated separation of the 
blister's wall from its ionising source (Baldwin et al. 1991; Schiffer 
& Mathis, 1974). At the edge of the blister the photon flux would 
be only ~3xl0^phot/s; as shown in Table 1, X-ray fluxes of order 
lO^phot/s produce a very cold X-region and as a consequence un- 
derestimates the BJ temperatures. 



5 CONCLUSIONS 

In this letter we have explored the possibility that X-ray ionised 
dense clumps embedded in lower density EUV-ionised gas may 
explain the discrepancy between ionic abundance determinations 
based on the analysis of RLs and CELs in H II regions. While 
the model can, in principle, simultaneously reproduce the observed 
ADF and BJ temperature in the Orion Nebula, the X-ray flux re- 
quirements do not match the observed X-ray emission. 

Independent from our current findings in the Orion Nebula, 
future observational studies may be used to empirically probe the 
X-E model by searching for a correlation between local ADF vari- 
ations with (measured or expected) local X-ray fields. If notwith- 
standing our pessimistic analysis above, this model were found to 
hold, the main implications for emission line abundance determi- 
nations would be that both CEL- and RL-determined abundances 
would require adjustments. CEL corrections are readily available 
using the simple method derived in this work. Furthermore, unlike 
the t" paradigm, the X-E model predicts that the error on the abun- 
dances derived from CELs is always smaller and easier to correct 
for than that on RL-determined abundances. This stems from the 
fact that RLs include contributions from both the X and E regions 
which have very different temperatures and ionisation structures. 
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The CELs on the other hands are mainly produced in the E re- 
gion and error on those is only introduced by the fact that the H^S 
line which is used in the abundance analysis contains contributions 
from the two regions. However this is easy to correct for by follow- 
ing the procedure outlined in this paper. 

In this letter we have mainly concentrated on the H II regions, 
however, larger ADFs have been recorded in some PNe (e.g. Wes- 
son et al. 2003, 2008; Ercolano et al. 2004; Liu et al. 2006). In 
recent years the large difference in the magnitude of the ADFs for 
PNe and H II regions has led several authors to conclude that dif- 
ferent effects may be at play in these systems (e.g. Garcia-Rojas et 
al. 2007). It is beyond the scope of this paper to discuss the abun- 
dance discrepancy problem in PNe in detail, however we specu- 
late that also in these system cold, quasi-neutral X-ray irradiated 
dense gas may the culprit for the high ADFs reported. X-ray emis- 
sion has indeed been observed in a number of PNe (e.g. Chu et 
al. 1997; Guerrero et al, 2005; Montez & Kastner, 2009), and it is 
well known that many PNe contain a dense equatorial disk, opti- 
cally thick to EUV radiation (e.g. Lester & Dinerstein, 1984). High 
density clumps have also been directly imaged in many PNe, like 
the Helix Nebula (O'Dell et al. 2004) and are likely to be present 
in many other objects. Furthermore, Tsamis et al. (2004) found that 
for a given PN the CEL C/O and N/O ratios were very similar to the 
RL C/O and N/O ratios, respectively. This lends further support to 
chemically homogeneous models like the X-E model, where CELS 
and RLs come from regions with the same heavy element abun- 
dances. PNe with their more easily constrained geometry, may in- 
deed, be a better laboratory to test the feasibility of this scenario. 
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